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•  ABSTRACT 

The  starting  point  of  this  paper  Is  a  recent  extension  by  A.K.  Theophf 
lou  of  the  Hohenberg-Kohn>Sham  (HKS)  density  functional  theory  to  ensembles 
of  systems  consisting  of  the  M  loMest  eigenstates,  equally  weighted.  As  In 
the  Hies  theory  the  key  quantities  are  the  exchange  correlation  energy, 
£^c^n(r)],Md  potent1aI,Vj^”  (r;Cn(r')J).  The  present  paper  provides  expres 
slons  for  these  quantities,  valid  for  systems  of  slowly  varying  density. 
Even  for  such  systems,  however,  there  are  essential  non-local  effects. 
Nevertheless  both  E^^g  and  v^”  can  be  calculated  in  terms  of  quantities 
characteristic  of  appropriate  uniform  thermal  ensembles.  This  theory  Is 
the  analog  of  the  ground  state  local  density  approximation  and  allows  cal¬ 
culation  of  excited  state  energies  and  densities. 


Permanent  address:  Department  of  Physics,  University  of  California 

Santa  Barbara,  California  93106. 


ror  \ 

i 

□ 

□ 

r 

.bility  Codes 

a  'd|or 

1 

In  its  original  formulation  ’  '  general  density  functional  theory  was 

a  ground  state  theory  for  non-relativistic  interacting  electrons  in  an 

external  potential,  v(r).  An  extension  to  ensembles  at  a  finite  temperature, 

3  2) 

0,  was  soon  developed  '  Hare  recently  the  formal  theory  was  extended 

4l 

to  “equl-ensembles"  consisting  of  the  lowest  M  states,  equally  weighted 
Both  inany*body  ensembles  are  characterized  by  appropriate  exchange  corre- 

a  '  u 

latlon  functionals,  F  [n(r)],E  [n(r)].  In  terms  of  these  the  exchange- 

xc  xc 

correlation  potential  of  the  Kohn-Sham  (KS)  equations  can  be  determined. 


5F  ®Cn(r*)]/6n(r)  L  thermal  ensemble  (la) 

XC  0 

<SExgCn(r')]/6n(r)(M  equi-ensemble  (lb) 


0 

and  the  ensemble  average  densities,  n(r),  free  energies,  «  ,  and  average 
energies,  E^,  respectively,  can  be. calculated. 

The  local  density  approximation  (LDA)  of  F  ”  has  been  previously  dis- 

XC 

5l 

cussed  In  the  present  note  we  develop  a  quasi-local  approximation  for 

M 

E_,  closely  related  to  the  LDA  for  thermal  ensembles. 

XC 

As  shown  by  Theophllou^^  the  average  density  n(r)  of  the  lowest  M 
excited  states^^  uniquely  determines  the  external  potential  v(r)  and  hence, 
implicitly,  by  means  of  the  Schroedinger  equation,  all  eigenstates 
For  every  n(r)  and  M  one  can  then  define  the  functional 


F^Cn(r)]  3  Av(i|»^,(T+U)i|ijjj)  , 
where  the  symbol  Av  has  the  meaning 


(2) 


(3) 


3 


T  and  U  are  the  kinetic  and  interaction  energy  operators;  and  the 

(m  «  are  the  M  lowest  eigenstates  corresponding  to  the  potential 

v(r)  which  reproduces  the  average  density  n(r)^^, 

n(r)  s  Av  n|jj(r)  .  (4) 

Using  the  functional  one  can  define  the  energy  functional 

I 

E^(p)Cn'(r)]  s  I  v(r)n'(r)dr  +  F^Cn'(r)]  ,  (5) 

whose  unique  tninimum  is  attained  when  n'(r)  is  the  correct  n(r)  and  has 

the  value  s  Av  E„. 

m 

This  minimization  can  be  carried  out  by  solving  appropriate  Kohn> 
Sham  (KS)  equations^^.  We  first  define 

T^[n(r)]  3  Ay  Tj_Jn(r)]  ,  (6) 

where  T"  Is  the  kinetic  energy  of  an  egu1>ensenib1e  jjf  non-interacting 

electrons  in  the  appropriate  external  potential  Vj(r)  yielding  the  given 

8) 

average  density  n(r)  i,e., 

Tg[n(r)]  »  n(r)dr  (7) 

where  E^  represents  total  single  particle  energy.  Now  we  can  define 
the  exchange-correlation  energy  functional  as 

E^”tn(r)]  a  F”[n(r)]  -  (T^[nir)]  .  \  |  dr  dr')  (8) 

and  the  effective  potential  by 


4 


M 

where  v^^(r)  is  given  by  (lb).  The  appropriate  KS  equations  are 

{-  y  V*  +  vj^^(r)  -  e^}  ♦'J(r)  -  0  .  (10) 

The  density  of  the  m'th  KS  eigenstate  is 
n^(r)  -  Z  f7|4”(r)|2  ^ 

where  f*J*  (■  1  or  0)  describes  the  occupation  of  the  i^  single  particle 
state  in  the  m'th  N«particle  state.  The  average  density  n(r)  is  then 
given  by  Eq.  (4).  For  example,  if  there  is  no  degeneracy  and  M  ■  2, 

N*1 

"(■•)■,  £  l♦?(r)l' •  •  (12) 

Eqs.  (10)  and  (4)  naist  be  solved  self-cons-lstently  for  n(r)  and  v|:^^{r), 
using  Eqs.  (11),  (9)  and  (lb).  Here  it  is  assumed  that  the  dependence  of 
E^^,  occurring  in  Eq.  (lb),  on  the  density  n(r)  is  known. 

It  remains  to  find  approximations  for  E^  and  v^  ,  in  the  spirit  of 

Aw  Xw 

the  LDA,  i.e.,  valid  for  systems  of  slowly  varying  density. 

THERMODYNAMIC  CONSIDERATIONS 

For  M  »  1,  a  non-degenerate  ground  state,  a  very  simple  and  useful 
approximation  for  e]  has  been  the  so-called  local  density  approximation 

Aw 

(LDA) 


EL[n{r)] 


(13) 


» L  ’I*  I  »y*-> '-» •wr^Tvryrxviv’r  -giy.-y- 1-- 

5 


.'  r-rf  r-  i  ^  .  •'■ 


where  exchange  correlation  energy  per  unit  volume  of  a  uniform 

electron  gas  of  density  n.  This  approximation  is  strictly  valid  only  when 
n(r)  is  a  slowly  varying  function  of  r.  However,  in  practice,  it  was 
found  to  yield  good  results  even  when  this  condition  was  not  satisfied. 

We  shall  now  generalize  this  approximation  for  an  equi-ensemble  with  arbi¬ 
trary  M  ^  1 . 

In  the  spirit  of  the  LOA,  we  shall  consider  systems  of  slowly  varying 
density  n(r).  Such  syst^ems  necessarily  occupy  a  large  volume  and  (unless 
n(r)  0}  contain  a  large  number,  N,  of  particles.  We  consider  both  the 
ground  state  of  such  a  system  (M  «  1)  and  equi-ensembles  of  the  M  lowest 
eigenstates.  Formally  we  may  consider  families  of  density  distributions 

n(r;a)  5  f(r/a)  a  »  a-,  ,82 .  (14) 


I 

i 


i . 

u  • 


r‘ 


( 
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where  f  is  a  given  function  and  a  Is  a  length  scale  parameter  which  becomes 
sufficiently  large.  We  denote  the  average  excitation  energy  of  the  equi- 
ensemble  by 

ae”  5  Av(E^)  -  E^  .  .  (15) 

As  a  the  spacing  between  excited  states  approaches  zero.  If,  as 

M  1 

a  -*■  «,  the  degree  of  excitation  as  measured  by  AE  /E  ,  remains  fixed, 
then  clearly  M 

In  such  a  limit  the  systems  can  be  described  by  the  principles  of 

thermodynami cs .  Accordingly  the  differences  between  a  canonical,  and 
91 

equi-ensemble  '  with  the  same  n(r)  and  same  mean  energy  become  negligible. 
We  can  write 


6 


M 


(16) 


where  S  1s  the  entropy.  We  denote  by  6  the  temperature  of  the  canonical 
ensemble  equivalent  to  the  equi-ensemble  with  M  states,  both  with  the 
same  density  n(r). 

We  now  study  the  relationships  between  these  two  ensembles  and  their 
corresponding  non-interacting  KS  ensembles.  With  the  aid  of  the  two  appro- 
prfate  exchange-correlation  potentials,  v”  and  v  ° ,  Eqs.  (la),  (lb),  we 

Aw 

M  A 

construct  the  two  effective  potentials  Eq.  (9),  and  similarly 

Next,  using  these  effective  potentitals,  we  solve  the  appropriate  KS  equa- 
tions  for  and  e”  ,  respectively,  and  form  the  two  non-interacting 

ensembles  corresponding  to  (M  or  S,  n(r))  and  (6,n(r))  respectively.  These 
two  KS  ensembles  are  not  Identical,  even  In  the  thermodynamic  limit.  The 
situation  Is  presented  In  Tableau  1.  By  construction,  the  temperatures  of 
the  canonical  real  ensemble  and  corresponding  KS  ensemble  are  equal,  9. 
Similarly  the  entropies  of  the  equi-ensemble  and  corresponding  KS  ensemble 
are  equal,  S.  However,  the  relations  between  entropy  and  temperature  are 
different  for  Interacting  and  non-interacting  ensembles  of  the  same  density, 
n(r).  This  Is  exemplified  by  uniform  ensembles  for  which  we  have 

e 

S(0)  =  l^^lr^de'  (Interacting)  (17) 

0 

and 


35(9) 


cJe') 

-=7T —  d0'  (non- Interacting) 

/  9 

0 


(18) 


with  unequal  heat  capacities  C  and  C^.  In  general,  referring  to  the 
tableau,  S  and  9  are  related  hv 


S  »  S(0;[n(r)])  .  (19) 

where  S(9)  describes  the  interacting  ensemble.  Similarly  and  9^ 
are  given  by 

S5  -  S5(9;Cn(r)])  (20) 

and 

I 

S  -  Sg(9^;Cn(r)])  (21) 

where  the  function  S^(9')  refers  to  the  non-interacting  ensemble. 

Using  these  relations  and  Eq.  (8),  we  can  write 

0 

_  -  Tj*[n(r)]  -  dr  dr' ,  (22) 

where  9  is  given  in  terms  of  M  by  Eq.  (19)  (with  S  •  k  log  M),  while  9 
is  given  in  terms  of  M  by  Eq.  (21).  Equivalently  we  can  write 

E,5[n(r)]  .  E/^[n(r)J  *  T®[n(r)]  -  T*®tn(r)]  ,  (23) 

■-here  the  thermal  exchange  correlation  energy  is  given  by 

E,*[n(r)]  =  -  T®[r.(r)]  -  dr  dr'  .  (24) 

For  systems  of  slowly  varying  density,  (23)  can  be  rewritten  as 

0 

E,"[n(r)]  .  I  e/^(n(r))dr  t  J  (t«(n(r)) -t5*(d(r))dr  ,  (25) 

d 

where  is  the  exchange  correlation  energy  per  unit  volume  of  a  uniform 

fl ' 

electron  gas  of  density  n  and  i-emDerature  d,  and  (n)  is  the  kinetic 


8 


energy  per  unit  volume  of  a  uniform,  non-interacting  electron  gas 

of  density  n  and  temperature  0 ' .  The  relations  determining  9  and 

9  become,  for  slowly  varying  density, 
s 


k  log  M 


^  9 

a  (n(r)dr 


(n(r) )dr. 


(26) 


9  6 

where  a  (n)  and  Og®(n)  are,  respectively,  the  entropy  per  unit  volume 
of  a  uniform  interacting  and  non-interacting  electron  gas. 

I 

Note  that,  since  9  and  6  depend  not  only  on  M  but,  implicitly, 

s 

on  the  entire  density  distributions  n(r) through  Eos.  (26)  the  super¬ 
scripts  9  and  9^  appearing  in  (25)  are  highly  non-local  functions 

M 

of  n(r') .  This  must  be  remembered  when  v  is  evaluated  by  taking 

xc 

M 

the  functional  derivative  of  e"  given  by  Eq,  (25)  (cf.  Eq.  lb). 

Let  us  re-write  Eq.  (25)  as 


xc 


(n(r) ] 


Je® (n(r)dr 


(n(r) )dr 


(27) 


where 


e®(n)  =  e®^(n)  +  t®(n),  (28) 

the  total  energy  per  unit  volume,  except  for  the  classical  electro¬ 
static  energy.  Then 


,M 


-M 


6E“  [n(r’) ] 


v“  (r; [n(r')]) 

Xw 


xc 
6n  (r) 


59 


60. 


3n  ^"^^n=n(r)  ’  6n(r) 


6n(r) 


■s,l 


{n(r') )dr 


e®(n(r'))dr' 


where  the  subscript  1  denotes  diffuentiation  with  respect  to 
the  temperatiire  eurgument.  Thus 


To  evaluate  56/5n(r)  at  constant  M  (or  S)  we  use  eq. 
which  gives 

similaurly  for  68^5n(r).  Thus  (29)  becomes 


/“c<^^?tn(r*)]  -  -  t/(n)^ 


n»n(r) 


Aa^(n)\ 

\  Sn  /n»n(r) 


/4(n(r-))dr’  4/»(n)\  /t® 

/aj(n(r'))dr-  V ' 


(30) 

(26) 

(31) 

(n(r'))dr 
d: _ 

®(n(r'))dr' 


For  convenience  I  now  describe  the  entire  cycle  of  solving  the  KS 
equations  for  an  equi -ensemble  of  N  states  In  the  quasl-LOA  (valid  for 
systems  of  slowly  varying  density). 

Consider  a  system  of  N  electrons  In  a  given  external  potential,  v(r}. 

The  objective  Is  to  calculate  the  average  density,  n(r),  and  average 
energy,  E^,  of  the  lowest  M  eigenstates^^. 

1.  One  requires  the  following  thermodynamic  functions]^ ^  for  homogeneous 
Interacting  and  non- Interacting  electron  gases ,  of  the  density  n  and 
temperature  6*.  (The  subscript  s  denotes  non-interacting 

•  t  » 

and  the  subscript  I  differentiation  with  respect  to 
temperature.) 

0 '  9 ' 

a)  The  entropies  per  unit  volume  o  (n)  and  o^  (n) , 

9  * 

(Eqs.  (26)}  and  their  temperature  derivatives, o^  (n)  and 

o®  (n) . 
s,l 

b)  The  exchange-correlation  plus  kinetic  energy  per 

unit  volume  e®  (n)  (Eq.  (28),  and  its  temperature 
9  * 

derivative  e^  (n) . 

c)  The  kinetic  energy  per  unit  volume  of  a  non-interacting 
system,  t®'(n)  (Eq.  (25), ff),  and  its  temperature  derivative 

2.  Begin  with  an  initial  approximation  to  n(r) ,  Determine  the 
corresponding  interacting  temperature  0  and  non-interacting 
(KS)  temperature  0^  by  solving  respectively  the  implicit 
equations 


11 


I 


k  log  M  »  /a® (n(r) )dr;  k  log  M  ■  fa_® (n(r) )dr 


3.  Construct  the  effective  one  particle  potential 


*  jlr-r'l 


where  v"  is  given  by  Eq.  (32) . 
xc 

4.  Solve  the  KS  single-particle  equations 

(-  7^*  +  VeJf(r)  -  e^)  »  0  .  (10) 

5.  Construct  the  M  lowest  non- Interacting  N-partIcle  wave-functions 
i|».  _(ra»l,...M)  and  calculate  their  average  density  n'(r)®^.  (See 
Eq.  (11)  ff.) 

6.  If  n'(r)  a  n(r),  then  the  original  n(r)  was  self-consistent.  If  not, 
repeat  steps  2-5,  starting  with  a  different  Initial  density  until 
self-consistency  Is  achieved. 

7.  wOw  determine  the  average  energy,  of  the  equi-ensemble  as 

follows.  Let  E^  ^(msl,...M)  be  the  energies  of  the  M  lowest 
KS  states.  Then 

^  ®s,«.  -  j/ ”|r-r'T’  ^ 

+  J{e®(n(r))  -  t^® (n(r) ) }dr  ,(33) 

We  now  add  some  remarks  about  the  thermodynamic  properties 
of  uniform  electron  gases,  listed  in  1)  above. 


•  a  -  j 
••  fc  F 

•\v\V 


VV  /  ^ 
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A  'k  -  -  ^  W  •  .'’■ 
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We  call  a  temperature,  0',  “low"  when  k9'<<Ep  where  Ep  is  a 


mean  Fermi  energy  of  the  ground  state  Ep  =<kp/2>=  j  3ir  <n(r)  '  >, 

the  bracket  denoting  an  appropriate  average.  In  this  regime, 
the  temperature  dependence  of  all  thermodynaunic  quantities  are  deter¬ 
mined  by  the  low  temperature  parameters  Y(n)  and  y.in)  characterizing 

S  '»  s*' 

the  linear  specific  heat  per  unit  volume? 


C®'(n)  =  Y(n)9'  ;  C®'(n)  -  Y„(n)0’, 

9  9 


(34) 


Some  calculations  of  the  thennodynafflic  functions  of  an  Interacting 

uniform  electron  gas  over  various  ranges  of  n  and  8'  have  already  been  re- 

ported^^*^^.  It  is  generally  believed  that  for  9'  ■  0,  has  been 

most  accurately  determined  (with  a  precision  of  order  0*1Z}  by  Monte  Carlo 
13 

methods  .  It  is  hoped  that  similarly  accurate  results  will  soon  become 

a  •  a  I  a ' 

available  for  the  finite  temperature  quantities  a  (n),  e”  (n)  and  v  (n). 

0  *  0 ' 

Calculation  of  the  non-interacting  quantities  a”  (n)  and  t”  (n)  is,  of 
course,  elementary. 


CONCLUDING  REMARKS 

The  reader  may  be  puzzled  by  the  rather  intricate  interplay,  in  this 
paper,  between  equi-  and  canonical  ensembles  of  different  temperatures  and 
different  entropies.  Indeed,  in  principle  a  knowledge  of  the  physical  prop¬ 
erties  of  canonical  ensembles  alone  determines  the  densities  n|jj(r)  (averaged 
over  any  multiplets)  and  energies,  E^(r)  of  all  eigenstates  For  example, 
let  us  write  the  partition  function  as 


Z(9) 


n(E)  dE  , 


(35) 
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where 


n(E)  «  E  5(E-E„)  .  (36) 

m  ™ 

Then  clearly  n(E)  is  the  Inverse  Laplace  transform  of  Z(0)  and,  by  Eq.  (3€), 
determines  the  positions  (and  multiplicities)  of  all  eigenvalues  E^j. 
Similarly  for  the  densities  n|||(r).  However  such  a  procedure  has  two  serious 

drawbacks.  Even  for  M  ■  2  it  requires  a  knowledge  of  Z(d)  for  a1J[  6  and  a 
calculation  of  all  single  particle  and  e^.  Secondly,  to  obtain  n(E) 
from  Z(0)  requires  an  anlytlc  continuation  into  the  complex  0-plane  Since 
Z(0),  for  real  0,  can  be  only  approximately  known,  such  a  continuation  may 
give  entirely  misleading  results. 

Why  then  do  we  not  deal  exclusively  with  the  equi -ensemble,  but  express 
both  E^-[n(r)]  and  v^  [n(r)]  by  means  of  thermal  quantities?  The  reason  is 
that  both  of  these  functionals  are,  even  for  systems  of  slowly  varying  den¬ 
sities,  highly  non-local.  There  is  no  simple  LDA  for  them,  i.e. 


E^leCMr)]  •  I  e^^(n(ry)dr 


^  L4(-.)  1 


n»n(r) 


not  possible.  (37) 


M 

For,  for  a  given  M,  the  local  contribution  to  E"  at  r  depends  not  only  on 
M  and  n(r)  but  also  on  the  entire  density  distribution  n(r'),  which  deter¬ 
mines  how  the  total  entropy,  S  ■  k  log  M,  is  apportioned  between  different 
volume  elements,  dr'.  On  the  other  hand  thermal  quantities  can  be  expressed 
in  the  form  of  a  simple  LOA,  e.q. 


kji'^ 
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E®j.Cn(r)]  *  I  ejj.(n(r))  dr,  (38) 

Q 

Since,  for  a  given  temperature,  the  local  contribution  to  depend  only  on 
the  local  n(r}  and  not  on  the  density  at  other  points.  We  can  take  advantage 
of  the  convenient  LOA  form  (38)  by  noting  that  an  equi>ensemb1e  is  equivalent 
to  a  thermal  ensemble  witti  the  same  n(r)  and  a  temperature  which  depends  on 
both  M  and  the  entire  density  distribution  n(r).  This  is,  of  course,  true 
of  interacting  real  ensembles  anH  of  non-interactina  KS  ensembles.  The  non- 
M  M 

locality  of  E^j.  and  enters  through  the  temperatures  of  the 
appropriate  corresponding  thermal  ensembles. 

Another  possibly  puzzling  issue  is  the  following.  It  may 
SewHi  questionable  whether  the  quasi -LOA  of  the  present  paper, 
derived  with  the  aid  of  thermodynamic  arguments  pertaining  to 
bulk  ensembles  with  very  dense  energy  spectra,  is  applicable  to 
the  lowest  few  states  of  the  system,  say  H  ■  1,  E  or  3.  We 
shall  now  explain  that  this  justified  question  is  of  the  same 
nature  as  the  question  whether  the  LOA  for  the  ground  state  is 
applicable  to  small  systems  of  2,  3  or  4  electrons. 

The  ground  state  LOA  is  logically  justified  only  for  systems 
of  many  electrons,  N>>1  ,  with  slowly  varying  density,  n(r); 
for  the  physical  assumption  underlying  the  integral,  Eq.  (13), 

Cn(r)]  is  that,  local  1  y,  the  electrons  can  be  regarded  as  a 
uniform  electron  gas.  Nevertheless,  the  ground  state  LOA  yields 
quantitatively  useful  results  for  systems  with  as  few  as 
2,  3  or  4  electrons  -  even  II 


.V'.- 


V  - 


Similarly,  the  present  quasi-LDA  is  logically  valid  only 
for  systems  of  many  electrons  whose  density  varies  on  a  large 
enough  length  scale,  a  (see  Eq.  (14)  and  for  M>>1 .  For  such 
systems,  as  already  mentioned,  the  spacing  between  excited 
states  approaches  zero  as  a-»-**  .  Therefore,  for  large  enough  a, 
there  exists  a  value  of  H  which  simultaneously  satisfies  the 
following  two  conditions: 


2.  (4e“o)/e^<6  f39) 


(See  Eq.  (IS)),  where  6  is  arbitrarily  small. 

In  view  of  the  first  condition  ,  (39),  thermodynamic 
considerations  such  as  the  equivalence  of  equi -ensembles  and 
canonical  ensembles  hold  for  M>N^.  On  the  other  hand,  because 
of  the  second  condition,  (39)  the  qua^i-LOA  yields 

for  M<M„ 

-  0 


with  arbitrarily  small  error.  This  is  in  fact  the  physically 


correct  result  under  the  ser-ond  condition  (39). 


How  useful  the  quasi -LOA  is  for  excited  states  when  the 
conditions  N>>1,  and  n(r)  slowyly  varying,  are  not  well 
satisfied,  remains  to  be  seen. 

By  successive  calculations  for  Increasing  M,  starting  with 
H  •  1,  the  excited  state  energies  and  densities  n^(r) 
(averaged  over  oiultlplets)  can  be  obtained  In  the  quasi -LDA  of 
the  present  paper. 
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